Atmospheric neutrino induced muons in the MACRO detector 



F. Ronga (for the MACRO collaboration) 

INFN Laboratori Nazionali di Frascati, 
P.O. Box 13 1-00044 Frascati Italy 

A measurement of the flux of neutrino-induced muons using the MACRO detector is presented. Different 
event topologies, corresponding to different neutrino parent energies can be detected. The upward throughgoing 
muon sample is the larger event sample. For this sample, produced by neutrinos having an average energy 
around 100 GeV, the ratio of the number of observed to expected events integrated over all zenith angles is 
0.74 ± 0.036stat ± 0.046s ys t ± O.lStheor- We have investigated whether the observed number of events and the 
shape of the zenith distribution can be explained by an hypothesis of — > v r oscillation. The best probability 
(17%) is obtained for sin 2 26 ~ 1.0 and Am 2 of a few times 10 -3 eV 2 , while the probability for the no oscillation 
hypothesis is 0.1 %. The other samples are due to the internally produced events and to upward-going stopping 
muons; the average parent neutrino energy is of the order of 4 GeV. The low energy data sets show a deficit of 
observed events similar to the one predicted by the oscillation model with maximum mixing suggested from the 
upward throughgoing muon sample. 



1. Introduction 

The interest in precise measurements of the flux 
of neutrinos produced in cosmic ray cascades in 
the atmosphere has been growing over the last 
years due to the anomaly in the ratio of con- 
tained muon neutrino to electron neutrino inter- 
actions. The observations of Kamiokande, 1MB 
and Soudan 2 are now confirmed by those of 
SuperKamiokande with larger statistics and the 
anomaly finds explanation in the scenario of 
oscillations (l). 

The effects of neutrino oscillations have to ap- 
pear also in higher energy ranges. The flux of 
muon neutrinos in the energy region from a few 
GeV up to a few TeV can be inferred from mea- 
surements of upward throughgoing muons |2j . As 
a consequence of oscillations, the flux of upward 
throughgoing muons should be affected both in 
the absolute number of events and in the shape 
of the zenith angle distribution, with relatively 
fewer observed events near the vertical than near 
the horizontal due to the longer path length of 
neutrinos from production to observation. 

Here the measurement about the high energy 
muon neutrino flux is presented, together with 
the first results on low-energy neutrino events in 



MACRO. [] 
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Figure 1. Sketch of different event topologies in- 
duced by neutrino interactions in or around MACRO 
(see text). In the figure, the stars represent the scin- 
tillator hits. The time of flight of the particle can be 
measured only for the Internal Up and Up Through 
events. 



2. MACRO as a neutrino detector 

The MACRO detector provides an excellent 
tool for the study of upgoing muons. Its large 
area, fine tracking granularity, symmetric elec- 
tronics with respect to upgoing versus downgo- 
ing muons and fully-automated analysis permit 
detailed studies of the detector acceptance and 
possible sources of backgrounds to the upgoing 
muon measurement. In addition, the overbur- 
den of the Gran Sasso Laboratory is significantly 
larger than that surrounding other experiments 
(Baksan and 1MB), hence providing additional 
shielding against possible sources of background 
induced by down-going muons. In our first mea- 
surement of upgoing muons |2| , we reported on 
a deficit in the total number of observed upgo- 



Figure 2. Distribution of the parent neutrino en- 
ergy giving rise to the three different topologies of 
events (see Fig. |l]), computed by Monte Carlo using 
the same cuts applied to the data. The distributions 
are normalized to one year of data taking. The aver- 
age energies of the three samples are about 4 GeV, 4 
GeV and 100 GeV, respectively. 



ing muons with respect to the expectation and 
also on an anomalous zenith angle distribution. 
In particular, too few muons were observed near 
the nadir. Here, we report on a much larger data 
set H which retains the same basic features as 
reported previously but with larger statistics. 

The MACRO detector is located in Hall B of 
the Gran Sasso Laboratory, with a minimum rock 
overburden of 3150 hg/cm 2 . It is a large rectan- 
gular box, 76.6 m x 12 m x 9.3 m, divided lon- 
gitudinally in six similar supermodules and ver- 
tically in a lower part (4.8 m high) and an upper 
part (4.5 m high). The active detection elements 
are planes of streamer tubes for tracking and of 
liquid scintillation counters for fast timing. The 
lower half of the detector is filled with trays of 
crushed rock absorbers alternating with streamer 
tube planes, while the upper part is open and con- 
tains the electronics racks and work areas. There 
are 10 horizontal planes in the bottom half of the 
detector, and 4 planes on the top, made of wires 
and 27° stereo strip readouts. Six vertical planes 
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Figure 3. Distribution of 1//3 for the full detector 
data set. A clear peak of upward muons is evident 
centered at 1//3 = —1. The widths of the distribu- 
tions for upgoing and downgoing muons are consis- 
tent. The shaded part of the distribution is for the 
subset of events where three scintillator layers were 
hit. 



of streamer tracking cover each side of the detec- 
tor. 

The scintillator system consists of three layers 
of horizontal boxes, with vertical boxes along the 
sides of the detector. The time (position) resolu- 
tion for muons in a scintillator box in this analysis 
is about 500 ps (~ 11 cm). 

Figure [l] shows a schematic plot of the three 
different topologies of neutrino events analyzed 
up to now: Up Through events, Internal Up 
events and Internal Down together with Up Stop 
events. Figure ^ shows the parent neutrino en- 
ergy distribution for the three event topologies. 
The requirement of a reconstructed track selects 
events having a muon. 

The Up Through tracks come from inter- 
actions in the rock below MACRO. The muon 
crosses the whole detector (E^ > 1 GeV). The 
time information provided by scintillator counters 
permits one to know the flight direction (time-of- 
flight method) . Almost 50% of the tracks inter- 
cept 3 scintillator counters. The average neutrino 
energy for this kind of events is around 100 GeV. 
The data have been collected in three periods, 
with different detector configurations. In the first 



two periods (March 1989 - November 1991, De- 
cember 1992 - June 1993) only the lower parts of 
MACRO were working. In the last period (April 
1994 - November 1997) also the upper part of 
MACRO was in acquisition. 

The Internal Up events come from v inter- 
actions inside the apparatus. Since two scin- 
tillator layers are intercepted, the time-of-flight 
method is applied to identify the upward going 
events. The average neutrino energy for this kind 
of events is around 4 GeV. If the atmospheric neu- 
trino anomalies are the results of oscillations 
with maximum mixing and Am 2 between 10 -3 
and 10~ 2 eV 2 it is expected a reduction in the 
flux of this kind of events of about a factor of two, 
without any distortion in the shape of the angu- 
lar distribution. Only the data collected with the 
full MACRO (live-time around 3 years) have been 
used in this analysis. 

The Up Stop and the Internal Down events 
are due to external interactions with upward- 
going tracks stopping in the detector (Up Stop) 
and to neutrino induced downgoing tracks with 
vertex in lower part of MACRO (Internal 
Down). These events are identified by means of 
topological criteria. The lack of time information 
prevents distinguishing the two sub samples. The 
data set used for this analysis is the same used for 
the Internal Up search. An almost equal num- 
ber of Up Stop and Internal Down is expected if 
neutrinos do not oscillate. The average neutrino 
energy for this kind of events is around 4 GeV. 
In case of oscillations we expect a reduction in 
the flux of the Up Stop events similar to the one 
expected for the Internal Up events, while we do 
not expect any reduction of the Internal Down 
events (having path lengths of the order of 20 
km). 

3. Upward throughgoing muons (Up 

Through) 

The direction that muons travel through 
MACRO is determined by the time-of-flight be- 
tween two different layers of scintillator counters. 
The measured muon velocity is calculated with 
the convention that muons going down through 
the detector are expected to have 1//3 near +1 



while muons going up through the detector are 
expected to have 1/(3 near -1. 

Several cuts are imposed to remove back- 
grounds caused by radioactivity in near coinci- 
dence with muons and showering events which 
may result in bad time reconstruction. The most 
important cut requires that the position of a 
muon hit in each scintillator as determined from 
the timing within the scintillator counter agrees 
within ±70 cm with the position indicated by the 
streamer tube track. 

When a muon hits 3 scintillator layers, there 
is redundancy in the time measurement and 1/(3 
is calculated from a linear fit of the times as a 
function of the pathlength. Tracks with a poor 
fit are rejected. Other minor cuts are applied for 
the tracks with only two layers of scintillator hit. 

It has been observed that downgoing muons 
which pass near or through MACRO may produce 
low-energy, upgoing particles. These could ap- 
pear to be neutrino-induced upward throughgo- 
ing muons if the down-going muon misses the de- 
tector (U . In order to reduce this background, we 
impose a cut requiring that each upgoing muon 
must cross at least 200 g/cm 2 of material in the 
bottom half of the detector. Finally, a large num- 
ber of nearly horizontal (cos0 > —0.1), but up- 
going muons have been observed coming from az- 
imuth angles (in local coordinates) from -30° to 
120°. This direction contains a cliff in the moun- 
tain where the overburden is insufficient to re- 
move nearly horizontal, downgoing muons which 
have scattered in the mountain and appear as up- 
going. We exclude this region from both our ob- 
servation and Monte-Carlo calculation of the up- 
going events. 

Figure |^ shows the 1/(3 distribution for the up- 
throughgoing data from the full detector running. 
A clear peak of upgoing muons is evident centered 
on 1/(3 = — 1. There are 398 events in the range 
— 1.25 < 1/(3 < —0.75 which we define as upgoing 
muons for this data set. We combine these data 
with the previously published data || (with 4 ad- 
ditional events due to an updated analysis) for a 
total of 479 upgoing events. Based on events out- 
side the upgoing muon peak, we estimate there 
are 9 ± 5 background events in the total data set. 
In addition to these events, we estimate that there 




Figure 4. Zenith distribution of flux of upward 
throughgoing muons with energy greater than 1 GeV 
for data and Monte Carlo for the combined MACRO 
data. The solid curve shows the expectation for no 
oscillations and the shaded region shows the 17% un- 
certainty in the expectation. The dashed line shows 
the prediction for an oscillated flux with sin 2 26 = 1 
and Am 2 = 0.0025 eV 2 . 



are 8 ±3 events which result from upgoing charged 
particles produced by downgoing muons in the 
rock near MACRO. Finally, it is estimated that 
1 1 ± 4 events are the result of interactions of neu- 
trinos in the very bottom layer of MACRO scin- 
tillators. Hence, removing the backgrounds, the 
observed number of upgoing throughgoing muons 
integrated over all zenith angles is 451. 

In the upgoing muon simulation we have used 
the neutrino flux computed by the Bartol group 
||. The cross-sections for the neutrino interac- 
tions have been calculated using the Morfin and 
Tung parton distributions set Si 0. These par- 
ton distributions were chosen based on the good 
agreement of the resulting gt compared to the 
world average at E v = 100 GeV. The propaga- 
tion of muons to the detector has been done using 
the energy loss calculation by Lohmann et al. || 
for standard rock. The total systematic uncer- 
tainty on the expected flux of muons adding the 
errors from neutrino flux, cross-section and muon 
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Figure 5. A) Probability contours for oscillation pa- 
rameters for Vy, — > ^ T oscillations based on the com- 
bined probabilities of zenith shape and number of 
events tests. The best probability in the physical regn 
is 17% and iso-probability contours are shown for 10% 
and 1% of this value (i.e. 1.7% and 0.17%). B) Con- 
fidence regions at the 90% and 99% levels calculated 
according to reference ||. Since the best probability 
is outside the physical region the confidence intervals 
regions are smaller than the one expected from the 
sensitivity of the experiment. 



propagation in quadrature is ±17%. This the- 
oretical error in the prediction is mainly a scale 
error that doesn't change the shape of the angular 
distribution. The number of events expected inte- 
grated over all zenith angles is 612, giving a ratio 
of the observed number of events to the expec- 
tation of 0.74 ±0.036(stat) ±0.046(systematic) 
±0.1 3 (theoretical) . 

Figure ^ shows the zenith angle distribution of 
the measured flux of upgoing muons with energy 
greater than 1 GeV for all MACRO data com- 
pared to the Monte Carlo expectation for no os- 
cillations and with a — » v T oscillated flux with 
sin 2 26 = 1 and Am 2 = 0.0025 eV 2 (dashed line). 

The shape of the angular distribution has been 
tested with the hypothesis of no oscillation ex- 
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Figure 6. Comparison between measured and ex- 
pected number of low energy events versus cos(6). 
The dashed line is obtained assuming neutrino os- 
cillation with the parameter suggested by the Up 
Through Sample. Note that in the second plot the 
flight direction is unknow 



eluding the last bin near the horizontal and nor- 
malizing data and predictions. The % 2 is 26.1, 
for 8 degrees of freedom (probability of 0.1% 
for a shape at least this different from the ex- 
pectation). We have considered also oscillations 

— > v r . The best \ 2 m the physical region of the 
oscillations parameters is 15.8 for Am 2 around 
0.0025ey 2 and maximum mixing (the best x 2 ls 
outside the physical region for mixing > 1 ). 

To test oscillation hypothesis, we calculate the 
independent probability for obtaining the number 
of events observed and the angular distribution 
for various oscillation parameters. It is notable 
that the value of Am 2 suggested from the shape 
of the angular distribution is similar to the value 
necessary in order to obtain the observed reduc- 
tion in the total number of events in the hypothe- 
sis of maximum mixing. Figure A) shows prob- 
ability contours for oscillation parameters using 
the combination of probability for the number of 
events and % 2 of the angular distribution. The 
maximum of the probability is 17%. The proba- 
bility for no oscillation is 0.1%. Figure |^B) shows 
the confidence regions at the 90% and 99% con- 



fidence levels based on application of the Monte 
Carlo prescription of reference ||. We plot also 
the sensitivity of the experiment. The sensitiv- 
ity is the 90% contour which would result from 
the preceding prescription if the data and Monte 
Carlo happened to be in perfect agreement at the 
best- fit point. The allowed regions are smaller 
than the one you could expect form the sensitiv- 
ity of the experiment. This is because the best 
probabil ity is outside the physical region. 

The same procedure applied for sterile neutrino 
|l(| gives 2% as maximum probability. 

4. The Low Energy Events 

The analysis of the Internal Up events is simi- 
lar to the analysis of the Up Through. The main 
difference is due to the requirement that the in- 
teraction vertex should be inside the apparatus. 
About 87% of events are estimated to be inter- 
actions. The preliminary uncertainty due to the 
acceptance and analysis cuts is 10%. After the 
background subtraction (3 events) 85 events are 
classified as Internal Up events 

The Internal Down and the Up Stop events 
are identified via topological constraints. The 
main requirement is the presence of a recon- 
structed track crossing the bottom scintillator 
layer. All the track hits must be at least 1 m from 
the detector's edges. The criteria used to verify 
that the event vertex (or stop point) is inside the 
detector are similar to those used for the Internal 
Up search. To reject ambiguous and/or wrongly 
tracked events which survived automated analy- 
sis cuts, real and simulated events were randomly 
merged and directly scanned with the MACRO 
Event Display. Three different events subsamples 
are considered according to the minimum num- 
ber of streamer tube hits. We present here the 
sample with at least 3 streamer hits (correspond- 
ing roughly to 100 gr cm~ 2 ). About 90% of the 
events are estimated to be i/ M CC interactions. 
The main background for this search are the low 
energy particles produced by donwn-going muons 
After background subtraction (5 events) 120 
events are classified as Internal Down and Up 
Stop events. 

The Montecarlo simulation for the low energy 



events uses the Bartol neutrino flux g and the 
neutrino low energy cross sections reported in || . 
The simulation is performed in a large volume of 
rock (170 kton) around the MACRO detector (5.3 
kton). The uncertainty on the expected muon 
flux is about 25%. The total number of events 
and the angular distributions are compared with 
the predictions in Table [l] and in Figure ^. The 
low energy samples show an uniform deficit of the 
measured number of events over the whole angu- 
lar distribution with respect to the predictions, 
while there is a good agreement with the predic- 
tions based on neutrino oscillations. 

5. Conclusions 

The upgoing throughgoing muon data set is in 
favor of — * v r oscillation with parameters sim- 
ilar to the those observed by Superkamiokande 
with a probability of 17% against the 0.1% for 
the no oscillation hypothesis. However the shape 
of the zenith distribution gives a maximum prob- 
ability of only 4.6%. This could be due to a sta- 
tistical fluctuation or to some hidden physics. We 
exclude effects due to the detector. 

This oscillation hypothesis is also consistent 
with the MACRO low energy data. A combined 
statistical analysis of the three data samples will 
be performed in the future when more statistics 
will be available. 
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Event Summary. The predictions with oscillations arc for maximum mixing and Am 2 = 0.0025eV r2 
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Predictions (Bartol neutrino flux) 


No Oscillations 
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